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1. Introduction 



Over the past few years the QCDSF Collaboration has focused on simulations of lattice QCD 
with 2 flavors of dynamical quarks. The real world consists, however, of Nf =2+1 light quarks (up, 
down and strange). We therefore extend our previous simulations to Nf =2+1 where we continue to 
investigate hadron and quark masses, weak matrix element, hadron form factors, moments of parton 
distributions as well as a variety of other key parameters of the Standard Model. Our ultimate goal 
is to bring the systematic uncertainties down to or below the experimental errors. 

JLQCD found an unexpected first-order phase transition in the strong coupling regime at rela- 
tively heavy quark masses when they employed the plaquette gauge action and the 0(a)-improved 
Wilson fermion action in three-flavor QCD simulations [Jp. Using an improved gauge action should 
give us significantly better control on the continuum extrapolation. Additionally, it is important to 
reduce somehow the chiral symmetry breaking arising from the Wilson fermion formulation. A 
well-known method to attenuate this symmetry breaking is adding a clover term. The UV filtering 
method, which involves replacing covariant derivatives in the fermion action by smeared descen- 
dents, is also becoming standard. We employ a tadpole improved Symanzik gauge action and 
stout link smeared Wilson fermions with a clover term. We also improve the algorithm to reduce 
simulation costs. 



2. The Action 



The tadpole-improved Symanzik action we use for Nf =2+1 simulations is 
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where the coefficients cq, c\ are taken from tadpole improved perturbation theory: 
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with Co + 8ci = 1, where uq = (|Tr(£/pi aqU ette)) 4 . We write ft = \cq. In the classical continuum 
limit uq — > 1 the coefficients assume the tree-level Symanzik values [|2p cq = 5/3, c\ = -1/12. 
We continue to use clover fermions with the action 



Sf = 2^{foMx)-KiKx)ul(x-m+yMMx-pL) 



(2.3) 
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but replace the gauge links U M in all terms of the fermion action except the clover term by stout 
links [g] 

U,^U tJ (x) = e iQ ^ x) U^x), (2.4) 

with 
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where V^(x) is the sum over all staples associated with the link. We take a - 0.1 and perform 1 
level of smearing, corresponding to a mild form of UV filtering [|J]. In this status report we present 
results where we used the tree-level value for the improvement coefficient, i.e. csw = 1, or used a 
value obtained from tadpole-improved perturbation theory: 

1 , , , 6 20m 2 . 

csw = -7 [1+^(0.0070628 1 + 1. 142004a -4. 194470a 2 )], g 2 = ^— . (2.6) 

Uq p 20uq - 8 

Note that in the future we will use the results presented in [Q]. 

This action has many advantages over our previously used one. In particular, due to UV filter- 
ing, it is expected to have better chiral properties [§] and smaller cut-off effects [Q]. One may also 
hope that the tadpole-improved perturbative value of csw is close to the non-perturbative value. 



3. The Algorithm 

The standard partition function for Nf =2+1 improved Wilson fermions is 

Z= \ DUDilfDipe~ s , 



(3.1) 

5 = S g (fi) + Sl(K[,Csw) + S s (K s ,C S w), 

where S g is a gluonic action, S / is an action for the degenerate u- and d- quarks and S s is an action 
for the strange quark. After integrating out fermions 

S = S g (fi) - ln[detM/M,][detMjM,F . (3.2) 

We first apply even-odd preconditioning: 

detM/M,cxdet(l + r^) 2 det!2je/, [detM f s M s p oc det(l + T s 00 )[detQlQ s p , (3.3) 

where 
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We then separate det Qj Q\ following Hasenbusch [| 
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Finally we modify the standard action to 



det Q ) Qi = det wj W t det — i, W=Q+p. (3.5) 
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where 



S l det = -2Tr log[l + r 00 (^)] , S s det = -Tr log[l + T 00 (k s )] , 
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(3.7) 
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We calculate S f r using the RHMC algorithm [Q] with optimized values for n and the number of 
fractions. We now split each term of the action into one ultraviolet and two infrared parts, 

S\jy = S g , Sir-i = S ! det + S S det + , SlR-2 = S l f 2 +Sf r . (3.8) 

In [ 1C ] we have introduced two different time scales [ |TI| ] for the ultraviolet and infrared parts of 
the action in the leap-frog integrator. Here we shall go a step further and put S\jy, Sir-i and Sir_2 
on three separate time scales, 

5t\ .. /6t) 



nr) = K- 2 (|) A- Vk- 2 (f)r. 



b = v^*\v q (JL-\vJ - 



2m\m,2i \m\m2) \lm\m2j ' 

where n r = t/(6t) and the Vs are evolution operators of the Hamiltonian. The length of the trajec- 
tory t is taken to be equal to one in our simulations. 



4. Test calculations 

We first tested our algorithm on small lattices of size 4 4 and 8 4 . Fi gure [I] shows the acceptance 
ratio and e~ AH for {3 = 7.2, tc; = k s = 0.1245 and csw = 1.0 for various simulation parameters. We 
discarded for themalisation the first 200 trajectories and then calculated for each choice of the 
simulation parameters about 1000 trajectories, e ~ AH should be equal to one within error and this is 
a good indicator of the correctness of the program. As seen in Fig. [|, we can keep high acceptance 
and e~ AH ~l by tuning parameters. Simulation results for the same parameters on a 16 3 32 lattice 
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Figure 1: The acceptance (left) and e A// (right) for various simulation parameters. K:(m\,rii2,p)- (3,3,0.1), 
(1,1,0.1), (2,2,0.1), (2,2,0.2), (2,2,0.3), (2,2,0.4), (2,2,0.5), (3,3,0.5). The other simulation parameters are 
fixed to n T =20, CG residual for Monte Carlo res mt =10~ 10 , CG residual for Molecular Dynamics res,„d=10~ 7 
and the rational approximation by n=2, 20 fractions and range [0.01,3]. 



are presented in table [jj The lattice spacing roughly corresponds to the so-called fine run of the 



4 



A status report of the QCDSF Nf—2+1 Project 



Yoshifumi Nakamura 



MILC collaboration [12] which is using the same gauge action. Figure ^ is a plot of k^, which is 
obtained by extrapolating the valence pion mass to zero, as a function of (romps) 2 . We see that k% 
for Nf -2+1 is lower than for Nf—2. Assuming that k% also has a rather mild ft dependence, we 
may consider this as an indication that our new action is much more continuum like. 



r /a 


amps 


amy 




4 


4.89(15) 


0.9089(20) 


0.9520(25) 


1.4828(53) 


0.134599(63) 



Table 1: ro/a, m P s, my, m;v and k% obtained from partially quenched calculations on a 16 3 32 lattice for 
P = 7.2, ki = k s = 0.1245 and c sw = 1.0. 
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Figure 2: k v c versus (r m PS ) 2 at j3 = 5.2 -5.4 for N/=2 and at p = 7.2 for Nf =2+1. 

Results for/3 = 7.2, ki - k s = 0.1335 with tadpole improved csw on a 16 3 32 lattice are presented 
in table ^[ They were calculated from 200 trajectories varying n of eq. (3.7), the number of fractions 
and the precision of the coefficients for the rational approximation. Our results confirm that double 
precision coefficients are needed to obtain the correct value of e~ AH . Note that the average plaquette 
value as well as the average minimum and maximum eigenvalue are consistent for the different 
choices of the algorithmic parameters. But since S y> may be O(10 8 ) (e.g. on large lattices) the 
algorithm may not be correct when using single precision coefficients. For the parameters shown 
in the last row of table ||, the ratios for the force contributions from the different terms in the action 
are 
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Ff r for n=4 is 40% smaller than for n=2. 



Ffr 

10, — ^ 10, 

FdeU 



Fde 



90. 



(4.1) 



5. Conclusion 



In this contribution we have presented the status of our Nf =2+1 project. We found indications 
for our action to be better than our previously used action. Furthermore, we tested the correctness 
of our algorithm. The performance of our program for matrix multiplication is about 20% of the 
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2.07(76) 
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Table 2: Simulation parameters and results for the value of the plaquette, the integrated autocorrelation 
time of the plaquette, e~ AH , acceptance, minimum and maximum eigenvalues of Q^Q on a 16 3 32 lattice for 
(3 = 7.2, ki = k s = 0.1335 and tadpole improved csw=L612. The parameters are the number n of eq. (3.7), 
the number of fractions (fr.) and the precision of the coefficients for the rational approximation (pr.). The 
other parameters are fixed to « T =60, mi=3, rri2=2>, res me =10~ 10 , res m£ /=10~ 8 and p=0.1. 

peak performance on the SGI Altix 4700. We are planning to implement better integration schemes 
and test other separations, for instance, 

S\jv = S g , SjR-i-S l fl , SiR- 2 = S f2 + S s fr + S s det + S l det . (5.1) 
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